Systematic Entomology

Phylogeny and taxonomy are related but not identical: phylogeny estimates common ancestry, whereas taxonomy must also provide diagnosable, stable, informative, and useful classifications. We examine this distinction using turtle and whale barnacles of the superfamily Coronuloidea, a group with an extensive fossil record, marked ecological specialization, a long history of morphology-based classification, and a recent topology-driven reclassification that merged recognized higher taxa to preserve monophyly. These changes reduced diagnosability and biological coherence, weakened historical continuity, and relied on placements insufficiently tested against morphological and phylogenetic evidence. Our reanalysis of published molecular datasets shows that critical relationships are sensitive to marker composition and taxon sampling. Together with morphological considerations on extant and fossil forms, these results indicate that topology alone is insufficient for coherent higher classification. We argue that higher taxa should be evaluated using complementary criteria: diagnosability, stability, utility, and monophyly, rather than by mechanically translating a preferred tree into Linnaean ranks.

Species-level delimitation within the genus Trichomycterus remains one of the main systematic challenges within the Trichomycterinae, particularly in lineages characterised by conservative external morphology, high apparent intraspecific variation, and historical diagnostic criteria based primarily on body proportions, colouration and a few meristic characters. In central Chile, Trichomycterus areolatus has traditionally been interpreted as a widely distributed and morphologically variable species, whilst T. maculatus, originally described from "Santiago du Chili", has remained subordinate to this broad conception without a modern phylogenetic reassessment. Here we reassess the specific boundaries of T. areolatus sensu lato using an integrative approach that combines complete mitogenomes, estimates of genetic divergence and comparative morphology of the cephalic laterosensory system associated with the neurocranium. Phylogenetic analyses reveal T. areolatus sensu lato to be non-monophyletic and identify a deeply divergent lineage, geographically coherent and attributable to T. maculatus. This lineage differs from restricted T. areolatus by extensive mitochondrial divergence, comparable to that observed between recognised species of Trichomycterinae, and by discrete characters of the cephalic lateral line system, primarily related to the continuity of the supraorbital canal and the arrangement of the associated pores. The congruence between mitogenomic, nuclear and morphological evidence supports the revalidation of Trichomycterus maculatus Valenciennes, 1846, and calls for a more restricted geographic circumscription of T. areolatus. These results demonstrate that the diversity of Trichomycterinae in central Chile has been underestimated, modify previous interpretation of the distribution of the species involved, and highlight the value of integrating mitogenomics and neurocranial/laterosensory characters into the taxonomy of morphologically conservative siluriform lineages.

Arcellidae is a family of testate amoebae within Arcellinida (Amoebozoa), comprising three recognized genera: Arcella, Galeripora, and Antarcella. Although species in the family have been studied for nearly two centuries, many historically described taxa and major morphological groups remain unsampled at the molecular level. Here, we provide a comprehensive review of Arcellidae and generate new cytochrome c oxidase subunit I (COI) sequences for arcellid species from Canadian peatlands, focusing on tall-shelled Arcella historically classified in section Altae sensu Deflandre. COI phylogenetic analyses recover a strongly supported monophyletic clade corresponding to North American representatives of Altae, providing the first molecular corroboration of this morphologically defined group. Within this clade, we redescribe Arcella leidyana based on modern material from Eeyou Istchee (Quebec). We further describe Galeripora purdoni sp. nov. from a calcareous fen in eastern Ontario, representing a novel terrestrial lineage within the genus, and redescribe Galeripora artocrea, which we transfer to Arcella based on congruent molecular and morphological evidence. Phylogenomic analyses of Arcellidae isolates from the Protist 10,000 Genomes Project reveal an additional deep lineage basal to Arcella and Galeripora. Together, these results highlight hidden diversity and demonstrate the importance of integrative approaches for resolving arcellid systematics and refining its classification.

Webspinners (Insecta: Embioptera) are an unusual insect order that are known for their subsocial behavior and prolific silk-production. Due to their unique foreleg silk glands, and spider-like ability to produce silk throughout their entire life cycle, webspinners are hypothesized to have evolved silk independently from other arthropod lineages. To date, there are no reference-quality genomes available for the order, preventing the study of their silk gene origination and diversification. Here, we assembled PacBio HiFi reference genomes and characterized the silk genes present in two webspinner species, Aposthonia ceylonica and Oligotoma nigra. The genomes reveal multiple full-length copies of the primary Embioptera silk gene, e-fibroin, that have undergone both ancestral and recent gene duplications within the group. For both species, all e-fibroin paralogs show the presence of complex repeat units consisting of multiple exons and introns that are remarkably homogenized across each gene. We also used CT-scanning of the internal silk glands to provide details concerning the localization of silk production in foreleg tarsi, and interspecific morphology. Article summaryThis study introduces the first high-quality genomes for webspinners, enabling new research on silk for evolutionary biologists and materials scientists alike. The authors sequenced two embiopteran species, Aposthonia ceylonica and Oligotoma nigra, to compare silk genes and gland structure using micro-computed tomography, an imaging method that shows internal anatomy in detail. They found multiple copies of the primary silk gene in both species that likely arose from multiple duplication events at different evolutionary times. These silk genes exhibit unusual gene structure with hierarchically organized repeat units that are highly homogenized within a gene. The findings show that silk genes have a complex evolutionary history in webspinners and provide a foundation for studying silk diversity within the order, and in the broader context of insect silk.

The floras of Australia and New Zealand are characterised by very high endemism of bryophytes. The highly diverse genus Lepidozia (fingerworts; family Lepidoziaceae) is represented by more than 20 mostly endemic species in these countries. The genus has been monographed in New Zealand, but issues surrounding species circumscription remain. This study focused on the Lepidozia ulothrix clade, with members shared between Australia and New Zealand. We performed phylogenetic analyses of a multilocus data set and applied several methods for molecular species delimitation. Our results support the conspecificity of L. serrulata and L. ulothrix, which is compatible with patterns of morphological variation observed in Tasmanian and New Zealand populations. The distribution of L. septemfida extends to the Wet Tropics Bioregion of north-east Queensland, where the species is represented by a distinct subspecies rather than a completely different species. The morphological differences between the Queensland and south-east Australian populations, as well as the inconsistency among the species-delimitation results, support the division of the species into two subspecies. Tasmanian and Victorian specimens previously attributed to the New Zealand endemic L. hirta belong to a distinct and unnamed lineage. Our study shows that integrative taxonomic approaches can uncover overlooked diversity and refine existing species circumscriptions, even in relatively well-studied Australasian lineages such as Lepidozia.

The Asparagaceae is an economically and ecologically important angiosperm (monocot) clade, containing iconic garden plants such as asparagus, hyacinth and agave, and arid-adapted keystone species like Joshua tree and dragon blood tree. The approximately 2,500 species in the family are spread across seven subfamilies: Agavoideae, Aphyllanthoideae, Asparagoideae, Brodiaeoideae, Lomandroideae, Nolinoideae and Scilloideae. Relationships among the seven subfamilies are generally unresolved and the sources of gene tree-species tree discordance are largely unknown. Investigating the biology underlying such discordance requires a robust sampling of both taxa and molecular markers, due to the extremely short branches that form the backbone of the Asparagaceae phylogeny. Using targeted sequencing of 1,726 low-copy nuclear genes, and a nearly complete sampling of genera, we interrogate relationships among all seven Asparagaceae subfamilies and investigate the underlying sources of discordance across the phylogeny. Analyses using different phylogenomic approaches (multispecies coalescence vs. nuclear gene supermatrix vs. chloroplast gene supermatrix) each revealed extensive discordance across the Asparagaceae backbone, which may be associated with rapid lineage diversifications that occurred in Africa shortly following the origin of the family ~49 Ma. Specifically, we found that two branches, representing ancestors of Aphyllanthoideae and Nolinoideae, may be responsible for much of the backbone discordance. The placement of Aphyllanthoideae, represented by one extant taxon (Aphyllanthes monspeliensis), is controversial with discordant topologies across many studies. However, multispecies coalescence-based analyses from this study strongly support a sister relationship between Aphyllanthoideae and a larger Agavoideae + Scilloideae + Brodiaeoideae clade. According to gene tree quartet frequencies and ABBA-BABA tests, we cannot reject incomplete lineage sorting (ILS) as the main source of discordance for the placement of Aphyllanthoideae; our results leave open the possibility of an ancient reticulation involving an ancestor of Aphyllanthoideae and an extinct basal Asparagaceae lineage. Regarding Nolinoideae, we found that conflicting support for two possible bifurcating clades (Asparagoideae + Nolinoideae vs. Asparagoideae + Lomandroideae) was associated with skewed support for alternative gene tree topologies, suggesting that phenomena other than ILS may explain the observed discordance. A quartet-based ABBA-BABA analysis, used to test for signatures of introgression, revealed a putative ancient hybridization event between two divergent, ancestral lineages which may have contributed to the origin of the Nolinoideae clade in Africa ~45 Ma. In sum, this study elucidates the long-troublesome backbone of the Asparagaceae phylogeny and advances our understanding of the evolutionary and biogeographical origins of all extant subfamilies. More broadly, this study expands our understanding of the various sources of gene tree discordance that may arise from natural processes associated with rapid lineage radiations across the tree of life.

Catenulids are free-living flatworms, common in eutrophic freshwaters such as ponds, ditches, or peatbogs, with most of the diversity described from tropical regions to date. Although the majority of the species have been described from warmer climates, most molecular studies have been done on specimens from temperate zones in Europe. We addressed this gap by sampling for exotic species in localities available in a temperate climate. In this study, we investigated catenulid diversity in the greenhouses at the University of Warsaw Botanic Garden and recorded two species known only from tropical areas (Stenostomum paraguayense and Suomina evelinae) and one exotic species recorded previously from a greenhouse in Poland (Stenostomum corderoi). Additionally, in the latter species, we provide evidence for environmentally induced coloration of sensory pits, which has not been reported thus far. We placed the collected species on a phylogeny using barcoding of 18S, 28S, and COI genes and retrieved paraphyly of the family Catenulidae, with S. evelinae forming a sister group to the genus Paracatenula, and hence we propose a revision of its systematic position. In total, we recorded six species, including three with a wide cosmopolitan distribution (C. turgida, S. grande and S. tuberculosum), and provided sequences for five of them, three of which had no previous molecular records (S. paraguayense, S. evelinae and S. corderoi). Thus, we confirm that greenhouses represent an important source of exotic species for taxonomic work on microscopic invertebrates.

Premise of the study: Doryopteris J.Sm. is a genus of dry-adapted ferns found widely across the southern hemisphere and the tropics. In the Hawaiian Islands, the genus is represented by four endemic species and one putative hybrid. Two of the endemic species (D. decora and D. decipiens) are widespread across all major islands, while two (D. angelica and D. takeuchii) are federally and state-listed endangered single-island endemics.We aimed to resolve the evolutionary relationships and biogeographic history of Hawaiian Doryopteris. Methods: We used a target-capture approach and the GoFlag 408 flagellate land plant probe set to sequence more than 400 nuclear loci for 40 Hawaiian and 29 non-Hawaiian individuals in Doryopteris and the related genera Lytoneuron and Ormopteris. After processing the sequences into alignments and single nucleotide polymorphism sites, we performed phylogenetic inference and population-level principal component analyses. Key results: We infer that Hawaiian Doryopteris are monophyletic and are nested within the pantropical species D. concolor. Within Hawai'i, clades and clusters of individuals correspond to geographical, geological, and climatological factors, more so than species identifications, with a few exceptions. We posit that Hawaiian Doryopteris represents an incipient radiation at the very earliest stages of divergence, and that budding speciation has led to non-monophyly of several taxa. Conclusions Our inference using genome-wide sequencing data sheds light on the complex evolutionary history of Hawaiian Doryopteris. Overall, the phylogeny of Hawaiian Doryopteris shows strong links with geographical, geological, and climatological processes and little with morphological evolution.

Understanding the relative roles of vicariance and dispersal in shaping diversity and distribution patterns is central to historical biogeography. In this study, we investigate the historical biogeography of the ancient centipede genus Digitipes Attems, 1930 from South and Southeast Asia. First, we determined the phylogenetic position of the genus Digitipes within the Order Scolopendromorpha (n=414) by assembling primary and published sequences (n=30) for two mtDNA markers (COI, 16S) and one nuclear marker (28S) using Maximum Likelihood and Bayesian inference. We further used single-locus and multi-locus coalescent-based species delimitation methods to identify putative species within the genus Digitipes. We then used three fossil calibrations to estimate divergence times in a Bayesian framework, and used the resulting time-calibrated phylogeny for biogeographic analysis in a likelihood framework (BioGeoBEARS). The genus Digitipes was monophyletic with strong support, with SEA lineages nested within the Indian clade and sister to the D. barnabasi species complex from the Western Ghats. D. pruthii, the Eastern Ghats species, was nested with the Western Ghats species clade. A single-locus mPTP-based species-delimitation method suggested the presence of 24 putative species, far exceeding the number of morphologically described species, indicating an underestimation of species diversity. Divergence time estimates suggest that Digitipes began diversifying around 126 mya (100-159 mya), affirming its Gondwanan origin. Time-stratified ancestral area reconstruction suggested that early vicariance, followed by jump dispersal and range expansion, shaped the distribution of the genus Digitipes in South and Southeast Asia. There was one dispersal event from India to Southeast Asia, following a transient land connection between them, around 50 mya, supporting the Out-of-India hypothesis. Additionally, three jump dispersal events and five range expansions explained diversification within peninsular India. Particularly, D. pruthii originated from a jump dispersal event from the Central Western Ghats to the Eastern Ghats around 37 mya. Our results highlight the importance of an integrative taxonomic framework to delineate hidden diversity and to obtain robust species hypotheses for testing biogeographic hypotheses.

Polytypic species with large ranges may harbor unrecognized diversity because taxonomy ranks populations differing subtly in plumage as subspecies. The Ruddy Foliage-gleaner (Clibanornis rubiginosus) exemplifies this problem. It ranges from Mexico to Brazil, with 15 subspecies, and forms a non-monophyletic complex with two congeners, yet its songs had not been compared. We measured ten spectral and temporal variables on 104 recordings covering 14 of 15 subspecies. Bayesian linear mixed models showed three song groups: eight subspecies west of the Andes share a single-note song, whereas Amazonian and Guianan populations add a short introductory note and sing longer, lower-pitched songs. Within this group, watkinsorum sings the lowest-pitched and longest song and is phylogenetically closer to C. cinnamomeigula than to its Amazonian neighbors. The third group is C. cinnamomeigula alone, a white-eyed taxon in an otherwise dark-eyed group. Its high-pitched, vibrato song resembles none other in the genus. One-note and two-note songs differ in kind without intermediates, and every two-note taxon sequenced to date falls in one clade, so C. rubiginosus is paraphyletic. We recognize four species, C. rubiginosus sensu stricto, C. cinnamomeigula, C. watkinsorum, and C. obscurus. This raises Clibanornis from five species to eight and divides its only polytypic species.

Recent technological advancements have enabled the rapid generation of high-quality genomes across the tree of life, often resulting in multiple reference genomes for clades of phylogenetic interest. These reference genomes are often used to reconstruct phylogenetically informative loci from short-read data of newly-sequenced species. However, this approach can introduce reference bias where the reconstructed loci have erroneous similarities to those of the reference genome. Since reference bias can seriously affect downstream analyses, it is important to assess its presence in phylogenomic datasets. In this study, we propose PhyloRBT (phylogenetic reference bias test) to detect reference bias by reconstructing each locus multiple times using different reference genomes and then measuring the phylogenetic correlation between these reconstructions and the corresponding locus from the reference genomes. We applied PhyloRBT to hundreds of BUSCO loci reconstructed from short-read data of nine Eucalyptus species using 34 different reference genomes. Across the nine species, we found that more than a quarter of the reconstructed loci had significant evidence of reference bias. Excluding putatively biased loci from species tree inference resulted in a species tree topology that is more consistent with expectations from previous studies. In conclusion, PhyloRBT offers a straightforward way to detect reference bias in individual loci, and to selectively remove those biased loci from downstream analyses.

Kashmir voles (Hyperacrius) are poorly studied fossorial arvicolines of the southern Himalayan slopes, with two species recognized. We sequenced Cytb from voles collected in northern Pakistan and combined these data with existing GenBank sequences, yielding 194 tips spanning nearly all arvicoline genera. We analyzed five data sets with maximum likelihood and compared external measurements among species using principal component and linear discriminant analyses. We recovered five well-supported clades of Hyperacrius separated by ~12-14% p-distance, which we treat as distinct species. We recognize clade 1 as H. wynnei, and populations referred to H. w. traubi are genetically indistinguishable from H. w. wynnei. We recognize H. zygomaticus, new rank (clade 2), and tentatively allocate clade 5 to H. brachelix. Clades 3 and 4 likely represent undescribed species. Morphometric analyses of external measurements clearly distinguish H. wynnei but not clades 3 and 4 (overall classification accuracy 70.3%). Phylogenetic analyses including nuclear genes placed Hyperacrius as sister to a clade containing Microtini, Arvicola, Lagurini, Dinaromys, and Ellobiusini, whereas mitochondrial analyses weakly allied it with Ellobiusini. Its traditional placement in Clethrionomyini is thus unsupported, although its exact position remains unresolved. Alticola is nonmonophyletic, with A. macrotis sequences nested within Clethrionomys, likely reflecting mitochondrial capture.

Among Madagascar primates, the sportive lemurs (family Lepilemuridae) have seen their species diversity increase from eight in 2005 to 26 in 2009 mostly by applying the phylogenetic species concept to DNA barcode data. Despite the genus being speciose, only one case of sympatry is known from northern Madagascar, where two sportive lemur species described based on low mtDNA divergence, Lepilemur ankaranensis and Lepilemur milanoii, were found to co-occur at the center of their joint distribution range. Here, to clarify the taxonomy of these two species and examine their sympatry, we apply an integrative taxonomic framework to genomic and morphological data from 84 individuals of L. ankaranensis and L. milanoii, encompassing their entire distribution range and the forest of Analafiana, beyond their southernmost limit. Using clustering, multivariate, and isolation by distance analyses, we find no evidence of a sympatric zone and show that despite clear genetic differentiation between regions, the genomic and morphological diversity of the L. ankaranensis, L. milanoii-Analafiana group is clinal and explained by geographic distance. These results clarify that L. milanoii is a junior synonym of L. ankaranensis and that the Analafiana forest population belongs to L. ankaranensis, extending its distribution. It further implies that the 'sympatric' zone, the Andrafiamena forest, hosts conspecific individuals with slightly differentiated mtDNA backgrounds, rather than slightly differentiated sympatric species. Lastly, we re-evaluate the IUCN conservation metrics of L. ankaranensis, which continue to qualify as Endangered (EN) under the B1ab(i-v) criteria.

Oceanic islands are natural laboratories for investigating the processes that generate, maintain, and redistribute biodiversity through time. An ideal flowering plant system to investigate these processes is Grewia L. (Malvaceae, Grewioideae), comprising ca. 283 species, of which nearly 30% are endemic to Madagascar, where they display considerable morphological diversity. Using Angiosperms353 target-capture sequencing for 171 accessions (119 Grewia spp.), we reconstructed the first densely sampled phylogeny of Grewia and performed divergence-time estimation, diversification analyses, and ancestral-range reconstruction. Our results support a Malagasy origin of the crown group around 27 million years ago, followed by extensive in situ diversification of Grewia. Repeated dispersal events from Madagascar to continental Africa subsequently contributed to the establishment of new African lineages, whereas colonization of the western Indian Ocean islands was more recent and resulted in anagenetic dispersals. Africa acted as the principal source of dispersal toward Arabia and into tropical Asia, suggesting the presence of a Miocene Afro-Arabian migration corridor under warm and humid conditions. Subsequent eastward dispersal from Asia to Australia is consistent with the stepping stone pathway provided by the Sunda and Sahul regions during the Miocene. This complex dispersal history was likely aided by fleshy fruits adapted for endozoochory. Our new evolutionary framework provides a robust basis for future taxonomic revision of this morphologically and ecologically diverse genus, and for identifying which characters have contributed to its exceptional diversification.

Speciation can be a slow process for long-lived organisms with large populations due to persistent incomplete lineage sorting, while hybridisation can further obscure species boundaries through admixture. In eucalypts, there is tremendous morphological and environmental variation with almost a thousand described species, subspecies, and hybrid taxa falling into multiple subgenera and intermediate sections. Multilocus molecular studies have generated concatenated species trees supporting deeper divergence between the Corymbia + Blakella + Angophora clade with the rest of eucalypts. However, the proportion of gene trees that recover the concatenated tree topology, including the taxonomic relationships among individual subgenera and sections, remains underappreciated. In this study, we assessed the predictive power of current eucalypt taxonomy - represented as a concatenated tree - in explaining topological variation across gene trees. To do this, we first extracted 1,187 BUSCO loci from short-read data of 701 samples representing roughly 500 described species using CAPTUS. We then inferred individual gene trees and a concatenated species tree using IQ-TREE2. Finally, we calculated the concordance factors for every described subgenus and section on the concatenated tree, and developed a taxonomic consistency test that tracks the taxonomic placement of each sample across gene trees. We show that, even though most taxonomic groups are monophyletic on the concatenated tree with 100% bootstrap support, several groups exhibit very low concordance, indicating that most gene trees do not recover their relationships as represented by the concatenated tree. These results highlight important considerations for current eucalypt taxonomic classification and provide a phylogenomic inference pipeline that can be applied not only to future eucalypt studies, but also to other evolutionarily complex groups of taxa.

Motivation: The Interactive Tree of Life (iTOL) is widely used to display and annotate phylogenetic trees, but managing its format-sensitive annotation files impede reproducible high-throughput analyses. Among the maintained Python packages and versions evaluated, none combined template generation, taxonomic monophyly assessment and iTOL batch operations. Results: PyiTOL validates inputs, generates 31 iTOL template schemas (22 accepted by the live batch uploader), performs LCA-based monophyly classification with nested-monophyly detection, sampling-completeness states and polyphyletic subgroup decomposition, plus API upload and session replay. On a topology-constructed benchmark, all calls matched prespecified labels for 4,389 groups; on a 700-genome tree, binary mono/non-mono calls agreed with ETE4 for 409 genera; 17,294 GTDB R232 genera were processed in about 17 s. Availability and Implementation: PyiTOL 1.0.3 (Python [≥]3.10; Linux, macOS and Windows) is MIT-licensed at https://github.com/ZengZichao/PyiTOL and archived with test data at Zenodo (https://doi.org/10.5281/zenodo.22106806).

Biodiversity aggregators such as GBIF provide unprecedented access to global biodiversity data, yet their representativeness remains uneven across space and taxa. This study examined the spatial and taxonomic structure of global vascular plant data available on GBIF. Six filters were applied to the GBIF vascular plant dataset, resulting in the removal of 54% of all records. Together, the filters explained more than 90% of the identified spatial issues, with duplicate and missing coordinates accounting for most of the variation. A higher number of occurrence records was associated with a greater number of spatial issues. Record distributions became progressively more even at finer taxonomic levels, from orders to species. The time series of occurrences for species, genera, and families increased sharply after 1800 and continued to rise, with no apparent stabilisation. Of the 824 ecoregions covered, 73 accounted for 72% of all occurrence records. These ecoregions spanned all continents but were strongly concentrated in Europe, followed by North America and Oceania. The analyses reveal four key patterns: (1) data volume is positively associated with spatial issues; (2) a small number of taxa account for a large proportion of records, whereas many are represented by relatively few; (3) occurrence data aggregated by GBIF have increased continuously since 1800; and (4) record coverage remains highly uneven across the world's ecoregions. These results highlight the substantial contribution of biodiversity data aggregators to expanding access to biological information while demonstrating the persistent spatial and taxonomic biases that shape their contents. Such biases should be explicitly considered when assessing data completeness and quality and when using aggregated occurrence records to infer global biodiversity patterns.

The Himalayan-Tibetan-Hengduan (HTH) region is the richest biodiversity hotspot for high-elevation plants. However, owing to its remote and physically challenging topography as well as the trans-national nature of the region, many taxonomic problems in the area remain unresolved, particularly in the Himalaya. This, in turn, has impeded our understanding of the assembly of its extraordinary high-elevation flora. Here, we resolve phylogenetic relationships and delimit species in a distinctive clade of honeysuckles that is endemic to the HTH, the Isoxylosteum clade of Lonicera, using restriction-site associated DNA sequencing (RADseq) and morphological data. Five species complexes of Isoxylosteum have standardly been recognized. Three of these complexes are highly variable and have been divided into several varieties or species each. Phylogenetic, population structure, and morphological analyses of leaf and floral traits from samples collected across the range of the clade support the recognition of five species, including a species that has most often been recognized as a variety of L. rupicola (L. rupicola var. minuta). Instead, we find that it is sister to L. spinosa. This is surprising because the geographic range of L. minuta is contiguous with the other varieties of L. rupicola in the northern Hengduan region but widely separated from L. spinosa whose range lies mainly to the west of the Tibetan plateau. On close examination we find that several morphological and ecological traits also support a closer relation of L. minuta to L. spinosa. None of the other eight previously recognized varieties and species were supported. Floral traits showed high discriminatory power, correctly classifying 89% of samples to species. By comparison, leaf dimensions classified species with 59% accuracy. Our results identify diagnostic morphological apomorphies for each recognized species and major clade and provide a revised taxonomic framework for Isoxylosteum.

Deep-learning classifiers for the fungal internal transcribed spacer (ITS) report accuracies above 90% and are increasingly proposed for environmental metabarcoding. We benchmark two pretrained models, a convolutional network and a transformer sharing a training corpus of 5.23M sequences, against two established k-mer methods on 5,222 identical queries, one per genus, evaluated on both full-length ITS and the ITS2 subregion that dominates environmental sequencing. The design favours the classifiers: queries are drawn from the same public dataset they were trained on and stratified by whether a querys genus lies in their own label space, recovered from the distributed models, while the reference the k-mer methods consult mirrors that label space and excludes the queries themselves. Even so, on full-length ITS both classifiers are outperformed by both classical methods at every rank and in both strata: SINTAX recovers the correct family for 92.0% of seen-genus and 70.3% of novel-genus queries and best-hit alignment against a 56,327-sequence reference for 92.3% and 67.2%, against 77.9% and 57.7% for the transformer and 76.5% and 54.3% for the convolutional network. A hierarchical logistic regression on k-mer counts, fitted in ten minutes to 1.07% of the MycoAI training corpus, also exceeds both and places novel genera better than either search method, and refitted on ITS2 it recovers 89.2% of seen-genus families on that amplicon against 88.6% for best-hit alignment, so neither learned classification nor the amplicon is what fails. Restricting the identical records to ITS2 costs the k-mer methods 3.5 and 3.7 percentage points of seen-genus family accuracy but costs the classifiers 49.6 and 58.0, reducing them to 28.3% and 18.5%. An ablation identifies the cause. Grafting each querys unaltered ITS2 between the flanking regions of a donor record from a different phylum returns the donors family for 34.4% of queries against the querys own for 4.2% in the convolutional model, and 63.7% against 0.8% in the transformer, from a baseline of 0.1% where no donor sequence is present. The models therefore read taxonomy principally from the flanking regions rather than from the ITS2 barcode, which explains the collapse and predicts the same failure for any subregion amplicon. Compounding this, on ITS2 the classifiers output probability all but ceases to separate novel from known genera (AUROC 0.541 and 0.503, the latter at chance, against 0.866 for alignment identity and 0.785 for the SINTAX bootstrap), so the failure is not detectable from the models own output. We recommend that reported accuracies for such models specify the amplicon region of the evaluation, state the length distribution of the training corpus, and include a same-query classical baseline.

Tephritid fruit flies of the tribe Dacini rank among the most economically destructive agricultural pests globally, with several Bactrocera Macquart and Zeugodacus Hendel species causing severe losses to fruit and vegetable production across South and Southeast Asia. Bangladesh harbors five dacine species of primary agricultural significance: Bactrocera dorsalis (Hendel), B. carambolae Drew & Hancock, B. zonata (Saunders), B. correcta (Bezzi), and Zeugodacus cucurbitae (Coquillett). Here we present a comparative mitogenomic and molecular phylogenetic framework based on 21 unique mitogenome records retrieved from NCBI GenBank, representing 19 dacine ingroup taxa and two outgroups. Direct parsing of the GenBank sequences confirmed the canonical complement of 13 protein-coding genes (PCGs) and two rRNA genes in all 21 records. Among the 14 Bactrocera ingroup taxa, genome size ranged from 15,273 to 15,977 bp. Whole-genome AT content ranged from 66.6% in Bactrocera tsuneonis to 82.2% in Drosophila melanogaster; the five Bangladesh-relevant dacine pest species showed tightly clustered AT content of 72.9-73.6%. A concatenated alignment of 13,600 nucleotide sites (13 PCGs + 12S + 16S rRNA) was analyzed by maximum-likelihood inference under the GTR+FO model (IQ-TREE v1.6.11; 1,000 ultrafast bootstrap replicates). The ML tree recovers the B. dorsalis complex (UFBoot = 79-100) and places B. correcta and B. zonata as a maximally supported sister pair (UFBoot = 100). This study provides a sequence-verified mitogenomic reference framework for molecular identification and pest surveillance in Bangladesh, and should be interpreted as a curated comparative baseline rather than a population-genomic analysis, as no newly collected Bangladeshi specimens were sequenced.